Capillary interactions in Pickering emulsions 
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The effective capillary interaction potentials for small colloidal particles trapped at the surface of 
liquid droplets are calculated analytically. Pair potentials between capillary monopoles and dipoles, 
corresponding to particles floating on a droplet with a fixed center of mass and subjected to external 
forces and torques, respectively, exhibit a repulsion at large angular separations and an attraction at 
smaller separations, with the latter resembling the typical behavior for flat interfaces. This change of 
character is not observed for quadrupoles, corresponding to free particles on a mechanically isolated 
droplet. The analytical results for quadrupoles are compared with the numerical minimization of 
the surface free energy of the droplet in the presence of ellipsoidal particles. 

PACS numbers; 68.03.Cd, 83.80.Iz, 83.80.Hj 



I. INTRODUCTION 

If colloidal particles get trapped at fluid-fluid interfaces 
they interact effectively via deformations of the interface. 
These so-called capillary interactions can easily be tuned 
by changing external fields and they depend sensitively 
on the shape of the particles. This makes them good 
candidates for designing self- assembling systems and pro- 
vides a convenient experimental playground for studying 
basic issues of statistical mechanics in two dimensions in 
the presence of long-ranged interactions [l|, 0|- On the 
other hand, recent experiments show that capillary 
forces between elongated particles floating on spherical 
interfaces can have important consequences for stabiliz- 
ing so-called Pickering emulsions, which are formed by 
particle-covered droplets (e.g., oil) in a solvent (e.g., wa- 
ter). 

Whereas considerable theoretical progress has been 
made in understanding capillary interactions at flat in- 
terfaces basic issues such as the balance of forces 
acting on the interface and the influence of the incom- 
pressibility of the liquid enclosed by spherical interfaces 
have not yet been fully resolved. Curved interfaces pose 
the additional difiiculty 0, Q that in the presence of ex- 
ternal forces the condition of mechanical equilibrium de- 
mands to either fix the center of mass of the droplet by 
an external body force or to pin the droplet surface, e.g., 
to a solid plate. The experimentally relevant issues of the 
boundary conditions at the plate and of their influence 
on the pair potential between capillary monopoles have 
been studied previously [§, [13] ■ Here, we study higher 
capillary multipoles and derive the corresponding effec- 
tive pair potentials for pointlike particles in terms of an 
expansion in spherical harmonics. By using a flnite el- 
ement method [ll| we compare these analytical results 



with those obtained by the numerical minimization of the 
surface free energy of drops in the presence of ellipsoidal 
particles. 



II. PERTURBATION THEORY 

We consider a spherical, full, and spatially fixed droplet 
of radius i?o with colloidal particles trapped at its surface 
with surface tension 7. We model the effect of particles 
in terms of an external surface pressure field Yi{VL) (see 
below) parameterized by spherical coordinates £7 = (6*, (j)) 
on the unit sphere. The equilibrium shape of the droplet 
subjected to the pressure field 11 follows from minimiz- 
ing the corresponding free energy functional J^[{ti(r2)}] 
expressed in terms of the dimensionless radial displace- 
ment of the interface v{yi) = {r{yi) — Rq) / R^. In the limit 
of small interfacial gradients |Vaw| <C 1 one has 0, [§| 



-{y avf - V'' - {-Kin) + ^Ji)v 

(1) 

-d^ is the dimensionless angular 



where Vq := eede + 

gradient on the unit sphere [l2| . The first two terms in 
Eq. ([T]) represent the surface free energy, the third term 
corresponds to the work done by the dimensionless exter- 
nal pressure 7r(ri) = n(r2)i?o/7 in displacing the inter- 
face and the fourth term serves to implement the volume 
conservation. The value of the Lagrange multiplier — /i 
follows from imposing the volume constraint 



(2) 
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jfi/Ra can be identified with the shift in the inter- 
nal pressure of the droplet with respect to the Laplace 
pressure 2"f/Ro of an unperturbed, perfectly spherical 
droplet. The stationary condition 5F /6v — leads to 
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the Euler-Lagrange equation 

^{Wl + 2)v{n) = 7r{n) + f,. (3) 

We expand both the deformation v{fl) and the pressure 
field 7r(r2) in spherical harmonics Fim(rJ) so that Eq. ^ 
turns into an infinite set of algebraic equations (T3| : 

[l{l + 1) - 2]vi.n, = TTi^ + fi5io, (4) 

where Z = 0, 1, . . . and m = —I, . . . ,1, and the spherical 
multipoles are defined as Xim ■= J d^l' X{n')Yim{^') for 
X{il) = v{fl) or X{n) — 7r(J7). The volume constraint 
in Eq. ^ implies vqo = and thus /i = — ttoo, which 
means that the internal pressure shift counterbalances 
the external pressure. According to Eq. (|4]) the I — 1 
components of the deformation v are undefined. This is 
consistent with the fact that those components describe 
translations of the whole droplet without any change in 
shape which do not change the free energy. On the other 
hand alU = 1 components of the external pressure tt must 
cancel, reflecting the condition of balance of forces acting 
on the droplet. Indeed, the multipoles 7ri_i, ttio, tth are 
proportional to the Cartesian components fx, fy, fz of the 
total force / acting on the droplet. Hence the condition 
/ = is equivalent to Trim — with to = —1, 0, 1. 



III. GREEN'S FUNCTION AND BALANCE OF 
FORCES 

As mentioned above we consider a pressure field in the 
form of the superposition of N pointlike forces fi pointing 
into directions fi^: 

N 

7r(f]) =^<Z,^(l],r!,), (5) 

1=1 

where q, := f^/i-fRo) and Sin,n,) = 6(9 - 9i)6{dp - 
sin 9i is the Dirac delta distribution expressed in 
terms of spherical coordinates. The pressure field as a 
whole must obey the condition of balance of forces which 
can be expressed as J^i' TTi^im = for m = — 1, 0, 1. The 
corresponding interface deformation can be written as 
the superposition of single particle contributions: 

JV 

v{n)=Y,Q^Gi^:^^). (6) 

i 

where G{il, is Green's function describing the re- 
sponse of the interface in direction f2 to a point-force 
applied in direction fl' and which, according to the con- 
straints of constant liquid volume and fixed center of 
mass obeys the equation 

~ (V^ + 2)G{n, f^') = E E YilMYimin'), (7) 

;>2 m=-l 
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FIG. 1: Angular configuration of the pressure distributions 
symbolically represented as gray patches centered at the di- 
rections Q.\ = (0, 0) (north pole) and Q2 on the unit sphere. 
The angles (pi and 02 represent orientations of the local coor- 
dinate frames x'y'z' and x"y" z" associated with the pressure 
distributions (see main text). 



where the right hand side is a modified Dirac delta func- 
tion j(r2 — r2') with the / = and I = 1 components pro- 
jected out. From integrating both sides of Eq. ([T]) over 
the unit sphere it follows that the I = component of G 
vanishes, which reflects the condition of constant volume 
J dil G{fl, ri') = 0. The I = 1 component of G also van- 
ishes, which can be seen by multiplying both sides of Eq. 
([7]) by a radial vector Sr, integrating over the unit sphere 
and using the fact that can be expressed in terms 
of spherical harmonics with I = 1 only. This refiects the 
condition that the center of mass of the droplet is fixed in 
space, which can be formulated as / dfi erG{fl, il') = 0. 

According to the above considerations the pressure 
field due to a single point-force violates the balance of 
forces acting on the droplet. Therefore, if one wants 
to cast the solution still into in the form given by Eq. 

with = 1, there must be an additional pressure 
field counterbalancing the single point-force. Inserting 
v{n) = G{n,n') into Eq. ^ and using Eq. dTj, one ob- 
tains the corresponding pressure field tt in the form 

■K{Vl)^5{n,VL')+TlCM{^:^'). (8) 

with the effective pressure ircMi^, — ~(3/(47r)) cos^, 
where 9 is the angle between the directions VI and U,' . The 
latter pressure must emerge in order to cancel the / = 1 
component of the external pressure -Kexti^) = 5{VL,Vl'). 
Therefore ttcm corresponds to a body force fixing the 
center of mass {CM) of the droplet without contributing 
to the deformation field v. 

The expression in Eq. ^ for the deformation under 
the action of point-forces can be generalized to the case 
of a pressure field ■Wf.xti^) due to an arbitrary continuous 
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force distribution: 



(9) 



The total pressure tt acting on the interface again has 
an additional component ttcm accomplishing the fixing 
of the center of mass, i.e., tt = Tr^xt + t^cm- One can 
express the free energy in terms of Green's function by 
integrating by parts Eq. ([l]), by using Eq. ^ and the 
divergence theorem, and finally by inserting v{VL) as given 
by Eq. © 13: 



iRl 2 



(10) 



IV. CAPILLARY INTERACTIONS 



We study the pressure T^ext localized around two dif- 
ferent directions J7i and 0,2- To this end we introduce 
the following decomposition: 



(11) 



where Ri {R2) denotes the rotation transforming the 
original coordinate frame xyz into the coordinate frame 
x'y'z' {x"y"z"), referred to as Oi (O2), associated with 
the pressure distribution tti {112) (see Fig. [1]). We use 
the parameterization in terms of Euler angles, in which 
an arbitrary rotation R can be composed of a rotation 
around the z-axis by the angle a, around the (rotated) 
y-axis by the angle /3, and finally around the (rotated) z- 
axis by the angle 7. Under the rotation R{a, f3, 7) of the 
coordinate frame the coordinates transform according to 
R^^, as indicated in Eq. pT|) . 

The total free energy can be written as = Fi,seif + 
F2,seif + Ai^, where i^,,,ei/ = -[ff /Ut:^)]\^{Ro/ a,) + 
0(1) is the self-energy of particle i @, which does not 
depend on the relative position of the particles on the 
droplet but depends on the particle size a^. For a particle 
located at the north pole this is defined as the smallest 
distance for which ni{9, (j)) —Q for all 9 > Ui/Ro and (p E 
[0, 27r); this corresponds to the solid angle circumscribing 
the particle-liquid interface. The interaction free energy 
AF is given by the cross-terms in Eq. (ITUl) with tt from 
Eq. (HI]): 
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X 



^AF = - I dn I dn'TTi{R^^n)G{n,n')x 



Ti2{R2^^') ^- Jdn Jdn' Y,^umYLiRI^n)x 



3,n 



where 
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1 

l{l -t- 1) - 2 



for 1^0,1 
for l>2 



(13) 



are the spherical harmonics coefiicients of G(r2, £7'). 
Without loss of generality we can take fii = such that 
the reference frame Oi coincides with the original refer- 
ence frame xyz (see Fig.[T|), which implies that i?i = 1. 

Spherical harmonics transform under the representa- 
tion of the group of rotations according to Yim{R^^^) ~ 
T!m'=-iDL',miR)Yim'{^), whcrc L*^, „ is the Wigner 
Z?-matrix and reads (adopting the convention used in 
Ref. [l|) 



Di,„(7?)=i?i,.„(a,/3,7) 



(14) 

where d^, „(/?) is known as the Wigner (small) d-matrix. 
We use the parameterization in terms of the orientations 
(pi and 02 of the coordinate frames Oi and O2, respec- 
tively, relative to the great circle connecting the points 
111 = and Q2 on the unit sphere (see Fig. [T]). The 
rotation R2 is then parameterized by the triad of the Eu- 
ler angles {2n — (pi, 6,4)2) ■ By using the orthogonality 
/ d^lYim{^)Yi*^,{n.) = SiifSmm' of the spherical harmon- 
ics, and due to the identities Yi*^{U) = (-l)™Yi._„j(ri) 

and 2?i-,™*(a,/3,7) = (-ir-'"'^l™',-™(«, Z?, 7) we fi- 
nally obtain 

^ I I 

7^0 



A. Limit of small particles 



i(m<f>i+m' tj>2) 



(15) 



Having in mind colloidal particles of radius a as sources 
of the effective surface pressure tt, we would like to relate 
the spherical multipoles nim, I > 0, m = —I, ... ,1 to pre- 
scribed forces, torques, and higher capillary multipoles 
defined for the case of a flat interface. If a <C i?o the inter- 
face can be treated as to be locally flat in the close neigh- 
borhood Afi of angular extent O{a/Ro) around each par- 
ticle separately. Assuming for reasons of simplicity that 
AQ is centered at the z-axis and approximating it by a 
circular disk D{a) of radius a in the tangent plane we can 
use the asymptotic form of spherical harmonics fisj 



Y, 



where 



Im \ 



Aln = 



-> I' 



i\+r. 



|m|t 



> 2l + l{l + n)\ 1 
47r (/-n)!2"n!' 



n > 0, 



(16) 



(17) 



in order to obtain 



TTim= I dn' 7T{n')Yi^{n') > 



i\m\+mj^ 



l\rn\ 



Rr 



\m\ + l 



(18) 
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where we have introduced the dimensionless moduh 



d^x' an(z') fz 



D(a) 



and the phase 



1 

- arg , 

lD{a) a 



(fx' an{z') 



1 



(19) 



(20) 



where rt > 0, z' denotes a complex number, and n(a;') = 
77r(r2'(a;'))/i?o, with x'{Q,') being a projection onto the 
plane tangent to the reference sphere at Vti. From Eq. 
(ITOl) it follows that the moduli Qn are the usual cap- 
illary multipoles defined at a locally fiat interface (see 
Refs. [SijilGj)- Accordingly, due to the conditions of force 
and torque balance on a flat interface @], Qq represents 
a "capillary monopole", i.e., the total external force act- 
ing on the particle, Qi a "capillary dipole", i.e., the total 
external torque. In the case of a free particle, Q2 repre- 
sents the lowest non-vanishing multipole. Furthermore, 
we note that for each particle i one can always choose 
the orientation (hj of the coordinate frame O, such that 



the phase of the order n vanishes, i.e., 4)i^n = 0. We can 
write the interaction energy lS.Fnn' for a pair of capillary 
multipoles Q„ and Q'.^, of arbitrary orders n > and 

> as a sum of terms such that in Eq. ((T5]) the sum 
over I is taken under the constraint / > max{2,n, n'} 
with m and m! fixed to to = ±n and to' = ±n'; this 
implies AF = n' =0 ^^nn' ■ Accordingly, the inter- 
action free energy /S.Fnn' for a pair of capillary multi- 
poles Qn and Qn' scales with the droplet radius Rq as 
Ai^n«' 7a"+ia'"'+Vi?^+"'. Note, that in the case of 
two monopoles {n — Q and n' — 0) the interaction does 
not depend on i?o but only on the moduh Qo and Qq. 

The Wigner d-matrix in Eq. (llSp can be expressed in 
terms of the Jacobi polynomials pjf'''^ (cos 0) [l^ : 



,{6) 



{I + m')l{l - m'y. 



{I + my.{l - my. 



1/2 



P, 



I— in' 



(cos 61). (21) 



J 



By using certain properties of the Jacobi polynomials (see, e.g., Ref. [Ta|) we finally obtain 



E 



(2^ + 1) 

7aa' a/7?o,a'/i?o^o' (-2)«+«'+in!n'!7r ^ ^ , f/-f2)(/-l) 

(i + n'y 



{i-ny 



(— 1)" cos(ri(/>i -I- n'02) ( cos - 



X < 



+ cos(n0i — n'(f)2) l^cos - 
(-l)"cos(n0i)P;"(cos^), 



71 +n 



sm- 



n +Ti 



P, 



j-f{n' —n,n' -\-n) 



(cos 6*) 



(cos 6*) 



, n > 0, n' > 0, 

n > 0, n' = 0, 
n = 0, n' = 0, 



(22) 



where Pi and P" are the Legendre and the associated Legendre polynomials, respectively. 



B. Effective interaction potentials 

For a first discussion, we evaluate the interaction en- 
ergy in Eq. ([22l) for a — a' and n' = n. In the case 
of two identical pointlike monopoles, n ~ 7i' = 0, one 
obtains [13] 



APoo(^) 



70^ 



47r 



1 4 

— \ — cos ( 

2 3 



2 cos In sin - 
2 



(23) 



which gives back Green's function G{0) = 
-AFoo(^)/(7a2g2) as given in Ref. [3. _ We 
note that according to Eq. ([231) one has G{9 
0) — --(l/(27r)) In^, which renders the deformation 



v{r) — — (l/(27r)) ln(r/i?o) due to a pointlike force, 
known for a flat interface, where r = 9Rq ^ Rq is the 
arc length (see the dashed line in Fig. [1]) and Rq plays 
the role of the capillary length. 

In the case of pointlike dipoles, n ^ n' = 1, one can 
evaluate the series involving the Jacobi polynomials by 
using their generating function (see the Appendix): 



70^ 



rAFii 



a/_RQ^O 

<p2)f+{e) + 



8tt \RoJ 



(24) 



where (f>i and (f)2 are the orientations of the particles as 
indicated in Fig. [T] and chosen such that the phases (j)i^n 
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defined in Eq. ([20)) vanish for n = 1, and where 

/ (0) 4sin2 - hi ( sin - V - sin^ - + 2, 

+ ^ ^ sin2(0/2) 2^2/3 2 ' 

(25) 

2 



/_(6') :=4 ( cos ^ ) In ( sin ^ | + ^ cos" 



(26) 



We note that the dependences on the angular separation 
9 and on the orientations 4>i , (/)2 do not factorize. There- 
fore, in order to minimize the free energy, the particles 
adopt orientations which depend on 9. According to Eq. 
(I?7)) below one can distinguish three branches of the free 
energy (see the dotted green lines in Fig. ^ correspond- 
ing to three different minimal orientations, out of which 
the one with the lowest free energy is the equilibrium one. 
Assuming that the rotational relaxation of the particles 
is much faster than the translational one, an effective in- 
teraction potential can be obtained by minimizing the 
free energy with respect to and 02: 



1 



{AFn(0",</)i,02)} 



{01.02} a/fio^O 

-f+{&) + f-{S), for 9<9o, 

'f+{9)-f-{9), for 9o<9< 

f+{9)-f^{9), for 9>h. 



Svr 



t t 



(27) 



where 6*0 arccos[l-2 exp(-10/3)] p; 0.38 and 6*1 ~ 1.83 
are the zeros of f-{9Q) — and f+{9i) — 0, respectively; 
the arrows indicate the minimal orientational configura- 
tion (e.g., t t corresponds to (i^i = 7r/2,(/i2 — 7r/2) etc). 

In the case of pointlike quadrupoles, n — n' — 2, one 
obtains 



1 



3Q2 
647r 



( A ) cos(2(^i 



a/Rg- 
4 



sm^{9/2) 



(28) 



Distinct from the case of dipoles, the free energy has 
a minimum for 4>i + 4>2 = kir, k = 0,1,... (so that 
cos(20i -I- 202) — 1), independent of 9. This corre- 
sponds to a monotonic attraction (see the solid blue line 
in Fig.©. 



V. DISCUSSION 

We discuss the influence of the curvature of the inter- 
face on the interaction potentials by comparing the ex- 
pressions for the curved interfaces (Eqs. ([M)) - (P5)) . solid 
lines in Fig. [5]) with the well established results for the 
case of a flat interface (see, e.g., Ref. [3, dashed lines 
in Fig. [5]), with the corresponding spatial separation d 
between the particles given hy d = Ro9. A striking dif- 
ference is the non-monotonicity for curved interfaces of 




FIG. 2: Reseated effective interaction potentials (Eqs. (|23|) - 
(|28p . solid lines) for two capillary multipoles of the same or- 
der n — 0,1, 2, each shifted as a whole such that AF„„{6 = 
tt) — 0. The dashed lines correspond to the following ex- 
pressions: [ln(e) - ln(7r)]/(27r) for n = 0, [l/vr^ - l/e^]/(27r) 
for n = 1, and 3[l/7r* — l/S*]/(47r) for n — 2, representing 
for comparison the case of a fiat interface with the equivalent 
spatial separation d = Ro9 of the particles. As expected, for 
9-^0 each full curve approaches its corresponding dashed 
curve. The green dotted lines correspond to three metastable 
branches of the free energy for n = 1 (see Eq. (127^ 1 with two 
relevant intersection points (at di ~ 1.83, where f+{6i) = 0, 
and ^0 — 0.38, where f-{6o) — 0) indicated by filled squares. 
In each interval [0, ^0], [^0,^1], and [^i,7r] the full green curve 
corresponds to the lowest of the three branches and as such 
describes thermal equilibrium. 



the free energy for monopole-monopole and dipole-dipole 
interactions, which in both cases leads to short-ranged at- 
traction and long-ranged repulsion. These peculiarities 
can be traced back to the presence of an external body 
force necessary for fixing the center of mass. For n = 0, 1 
the expected reduction of the results for curved interfaces 
to those for a flat interface occurs only for very small an- 
gular separations. This is not the case for quadrupoles 
(n = 2), for which at all separations there is no qualita- 
tive discrepancy and only a minor quantitative difference 
from the flat limit. 

Next, we compare our approximate analytic theory 
with the numerical minimization of the free energy for 
floating ellipsoidal particles with semi-axes a and b, which 
we performed by using a finite element method [11,] . The 
numerical procedure consists of minimizing the free en- 
ergy functional 

J'[{r(ri)}, hi,h2,ipi,ip2;0, cj)i,cj)2,9p, a, b, Vi, A] = 

= "/Sig - 7 ^ cos9p,,Spu - HV - Vi), (29) 
1=1,2 

with respect to the shape of the interface {r{n)}, immer- 
sions of the particles hi,h2, and their tilts ipi,4'2 (rela- 
tive to planes perpendicular to the directions fli and fl2, 
respectively, see the inset in Fig. while the angular 
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FIG. 3: Rescaled effective interaction energy AF^f^ for two 
ellipsoidal particles of aspect ratio a/b = 3, contact angle 
9p — 2tt/3, various droplet radii Ro, and two different ori- 
entational configurations: tip-to-tip and side-to-side (as vi- 
sualized by the snapshots from the finite-element free en- 
ergy minimization). In view of large relative numerical er- 
rors for small values of AF^f^ the zero of the free energy is 
taken as the mean value of the data points with 6 > 1.5 
(the error bars are comparable or smaller than the size of 
the symbols). The thick solid line corresponds to Eq. (|28p 
with — cj>2 = or (^1 = (f>2 = 7r/2 (in either case Eq. 
(|28[l yields the same free energy). The capillary quadrupole 
of each of the particles is approximated by the expression 
Q2 = 2nAu/a; this relation is valid for a flat interface [oj 
at large interparticle separations, with the undulation Au of 
the interface taken to be the one at the surface of a single 
ellipsoidal particle. As an approximation for a droplet we 
adopt the same relation, Q2{Ro) ~ 2tiAu{Ro) / a, with the 
undulation Au = max^ r{9 = a/ Ro, (p) — min^ r{6 = a/Ro, (j)) 
of the interface at the bottom part of a particle centered 
at the north pole of the droplet (see the green part of the 
upper inset). This undulation has been calculated numeri- 
cally for three droplet radii: Au{a/Ro = 0.599)/a = 0.028, 
Au{a/Ro = 0.5)/o = 0.031, Au{a/Ro = 0.375)/a = 0.034. 
Rf) is the radius of the droplet without the particles. The 
immersion hi and the tilt angle tpi of the long axis of each 
particle i with respect to the plane perpendicular to the di- 
rection Q,i are subjected to minimization. The colored lines 
end at small values of 6 because the particle surfaces touch 
each other for nonzero 9. 



separation 6 and the tangential orientations (f>i , (f>2 are 
kept fixed; Sig is the area of the liquid-gas interface, and 
Spi.i and 0p^i are the area of the particle-liquid interface 
and the contact angle, respectively, of particle i. The last 
term in Eq. (1291) ensures conservation of the liquid vol- 
ume Vi with A as the corresponding Lagrange multiplier. 

The effective potentials (Fig. [3]) obtained numerically 
for prolate spheroids (i.e., elongated spheres) with aspect 
ratio a/b ~ 3 and droplet radii a/Ro = 0.599,0.5,0.375 
differ for particles oriented tip-to-tip {(pi — (j)2 = 7''/2) 
and side-to-side (01 = 02 — 0). This signals the 
importance of capillary multipoles higher than only a 
quadrupole, because the point-quadrupole approxima- 



tion yields the same free energy for both orientations 
(see Eq. (^5]) ). This discrepancy is not surprising be- 
cause the longer particle axes (which set the length scale 
a) are actually not very small compared with the droplet 
radius Ro and thus in this case the point-particle limit 
is a rather crude approximation. Our numerical analysis 
is restricted to such small droplets due to the limited ef- 
fectiveness of the finite-element method for larger drops 
caused by too large numerical errors. For the same rea- 
son we consider only particles with a large aspect ratio 
[a/h — 3) for which the deformation of the interface is 
sufficiently large so that the numerical noise does not 
smear out the angular dependence of the free energy. 
Nonetheless, even in this regime our analytic theory cap- 
tures the main qualitative features of the effective inter- 
action potential between fioating ellipsoids. Presumably 
the reason for this is that even for large particles the 
droplet shape deviates only slightly from a sphere. (Due 
to a fixed volume of the liquid larger particles cause a 
larger displaced liquid volume Vdispi ; primarily this leads 
to an increased effective radius ^ ^Vi -\- Vdispi of the 
droplet and thus generates only a constant in the de- 
formation field V.) Accordingly, away from the parti- 
cles our perturbation theory remains applicable. Thus 
even for large particles the deformation of the droplet 
can be approximated by the superposition of contribu- 
tions from pointlike quadrupoles substituting the actual 
particles. The additional dependence on the orientations 
revealed by our numerical calculations, which goes be- 
yond the quadrupole approximation, resembles the well 
established case of a fiat interface, which has been ana- 
lyzed experimentally flS^ and theoretically pj| in terms 
of a multipole expansion. Those studies showed that 
the point-quadrupole approximation, within which tip- 
to-tip and side-to-side configurations have the same free 
energy, is reliable only either for particles with a very 
small eccentricity (& ~ a) or at large spatial separations 
d a. Otherwise, one should rather consider an ex- 
pansion in terms of elliptic coordinates, from which it 
follows that the tip-to-tip configuration is energetically 
more favourable. This prediction agrees with the exper- 
imental observations 18]. Our calculations indicate that 
the same configuration is the preferred one also in the 
case of spherical interfaces. However, we have not ana- 
lyzed all the possible orientations. It is noteworthy that 
optimal orientations other then tip-to-tip and side-to-side 
have been observed at a flat interface for a pair of ellip- 
soidal particles of different sizes poj . 

We can also discuss the effective capillary potential 
F^li {a, 01 ; 6*0) for a single capillary quadrupole on a ses- 
sile, spherical caplike droplet forming a contact angle 
with the supporting substrate, where a is the angular po- 
sition of the particle (a < 0o) measured from the droplet 
symmetry axis and relative to the center of the sphere 
circumscribing the droplet, and 0i the orientation of the 
particle. As before the tilt angle tpi is taken to be the op- 
timal one. In the case Oq = tt/2 the method of images can 
be applied, analogous to the case of monopoles [9]. The 
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substrate acts like a mirror and the free energy of a single 
particle equals one half of the energy of two particles on 
a full droplet, i.e., 

AFW(a, </)i; 00 - vr/2) := F^li>{a, 0i; 0o = ^2) 
-F^li>iO,^,;9o = 7r/2) = 

= ±^AF22{TT-2a,(j>l,TT-(f>l), (30) 

where the + sign corresponds to a free contact line at the 
substrate and the — sign to a pinned one. (Since the par- 
ticle is force-free, there are no additional images, distinct 
from the case of monopoles, see Ref. Q.) For a free con- 
tact line Eq. ([5(7|) holds in general, i.e., for arbitrary par- 
ticles (not necessarily ellipsoidal) . In such a case AF22 is 
replaced by the actual interaction energy AF^'^^ between 
the particles (possibly containing also higher multipoles) . 
On the other hand, for a pinned contact line Eq. ([50)1 
holds only for pointlike particles in conjunction with Eq. 



For a pointlike quadrupole, according to the expression 
for AF22 in Eq. (^5)) . the dependence of AfI:^ on the ori- 
entation 01 drops out so that within this approximation 
the particle is always attracted to a free contact line and 
repelled from a pinned one. Thus the effective potential 
is given by the black line in Fig. |3] with 6 = w — 2a (with 
an additional minus sign in the case of a pinned contact 
line) . 

In the case of an actual ellipsoidal particle the free en- 
ergy is no longer given by the rhs of Eq. (^5)) so that the 
dependence on (fix does no longer drop out but becomes 
important in the vicinity of the contact line at the sub- 
strate. As stated above, if this contact line is free (still 
assuming = '''/2), due to symmetry arguments the free 
energy of a single particle on a sessile drop equals one half 
of the energy for two particles on a full droplet (as pre- 
dicted by the method of images for the more restrictive 
case of pointlike quadrupoles) , for which the full numer- 
ical results are plotted in Fig. [31 Accordingly, for a fixed 
angular position a = 7r/2 — 0/2 of the particle its orienta- 
tion perpendicular to the contact line (solid lines in Fig. 
[3]) is energetically more favorable. However, the global 
minimum of the free energy corresponds to the particle 
touching the substrate (i.e., a — > 7r/2 so that — > 0) and 
being oriented parallel to the contact line (dashed lines in 
Fig. |3]for < 0.5). However, we cannot judge if the lat- 
ter configuration remains the optimal one for arbitrary 
aspect ratios and contact angles at the particle. (This 
and analyzing the influence of the contact angle at the 
substrate actually call for a separate numerical study.) 

If the contact line at the substrate is pinned, the global 
free energy minimum corresponds to the particle being 
positioned at the drop apex. In such a case, due to a 
large distance from the contact line with the substrate, 
the pointlike quadrupole approximation remains reliable. 
One can estimate the depth of this effective confining po- 
tential as follows. According to Eqs. ((50)) and (1^51) with 



Q2 = 2'KAu/a and assuming Au/a ~ 0.05 (which ap- 
plies to an ellipsoidal particle of aspect ratio a/b = 3 
and contact angle 9p = 27r/3, see the caption of Fig. 
121), the displacement of the particle, e.g., from a — 
to a = 7r/4, leads to a free energy increase of approxi- 
mately ksTx {a/ fim)'^ for a droplet as large as i?o/a — 10 
(assuming a surface tension 7 ~ 0.07N/m which corre- 
sponds to a water droplet in air). Hence, one could ex- 
perimentally measure the capillary confining potential, 
e.g., by observing the Brownian motion of micron-sized 
particles on sessile drops. Finally, we emphasize that 
the effects of confinement should be also taken into ac- 
count in determining the effective pair potentials between 
particles on sessile drops (for the case of monopoles see 
Ref. |l3l). In particular, the confinement leads to closer 
packing of particles on a droplet, which should not be 
mistaken for an additional attractive pair-interaction. 
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VII. APPENDIX: CALCULATION OF SERIES 
CONTAINING JACOBI POLYNOMIALS 

In Eq. ([22I) for n = n' = 1 the following series appears: 



(j + l)(j+2)(2j-f3) . 



where {ai,pi) = (2,0) or (0,2), and for n = n' = 2: 



(31) 



St'''\x) ^(j + 2)(j + 3)(2j + 5)Pj''^^'-^\x) (32) 

j>0 

where {(J2,P2) = (4,0) or (0,4). Those series can be 
evaluated using the generating function [2l|: 



g^^'P^x^z) 



1 



R{1- z + Ry{l + z + R)f^ 



T.-^^^ (33) 



20-+P 



where R — \/l — 2xz + and \z\ < 1. One finds the 
relationships 



s['''-"'\x) ^ 2 fdyy^ [' dz 
Jo Jq 



d a _2 9 



X z ''^^z^''^-^z- — z 
oz oz oz 



2''^+P^g<-''^'P^\x,z) -1] (34) 
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and and 
which can be evaluated with the resuhs 

(36) 

.(0,2), - ^l-x\ 20 



S!f'^\x)^0. (39) 



5r^(2;) = -21n(^^j--, (37) 

S2*''''\x) = ^ , (38) Taking x = cos 6 and inserting these results into Eq. (|22|) 



(1 - ' finally yields Eqs. ([Mjl -iP 
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